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We show how generic off-diagonal cosmological solutions depending, in general, on all 
/— ^ , spacetime coordinates can be constructed in massive gravity using the anholonomic frame 

deformation method. Such metrics describe the late time acceleration due to effective 
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cosmological terms induced by nonlinear off-diagonal interactions and graviton mass and 
include matter, graviton mass and other effective sources modelling nonlinear gravita- 
tional and matter fields interactions with polarization of physical constants and defor- 
mations of metrics, which may explain certain dark energy and dark matter effects. 

Recently, a substantial progress in quantum gravity was made when de Rham 
- -, , and co-authors had shown how to eliminate the scalar mode and Hassan and Rosen 

Q^' established a complete proof for a class of bigravity / bimetric gravity theories.'^ 

5^ , The second metric describes an effective exotic matter related to massive gravitons 

W)| and does not suffer from ghost instability to all orders in a perturbation theory and 

away from the decoupling limit. 

We study modified massive gravity theories constructed on a pseudo- 
^ , Riemannian spacetime V with physical metric g = {g^ij/} and fiducial metric as 

^^3 ' we shall explain below. The action for our model is 

^'i . . . ^. (.^) 

^S , where /i = const is the mass of graviton, R is the scalar curvature for an " auxil- 

CO ' iary" (canonical) connection D uniquely determined by g via conditions^ 1) Dg = 0, 

and 2) its h- and w-torsions are zero for the nonholonomically induced torsion T 
completely determined by g for a conventional splitting N : TV =/iV vY. The 
"priority" of D is that it allows to decouple the field equations in various gravity 
^ , theories and construct exact solutions in very general forms. We can impose addi- 

tional nonholonomic (non-integrable constraints) when D.s^^^ — > V and R — > R, 
where R is the scalar curvature of the Levi-Civita (LC) connection, V, and it is 
possible to extract exact solutions in GR and various torsionless modifications. We 
use the units when h = c = 1 and the Planck mass Mpi is defined via Af|,j = I/SttG 
with 4-d Newton constant G. 

The equations of motion for modified massive gravity theory, MGT, ([T]) are 
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(9j^/)R^, - -/(i?)g,,. + /i^X^, = Afpf T^„ (2) 



where R^^ is the Einstein tensor for g^,^ and D and T^,^ is the standard matter 
energy-momentum tensor. For D — > V, we get Hfn, — > Rf^, with a standard 
Ricci tensor Rf^i, for V. The effective energy-momentum tensor X^i, is defined in a 
"sophisticate" form by the potential of graviton iY = Z^2 + «3^3 + ^4^4 , where a^ and 
0:4 are free parameters. The values U2 , W3 and U4 are certain polynomials on traces of 
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some other polynomials of a matrix IC^ ^ '^u ~ ( v 5 ^^ ) for ^ tensor determined 
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by four Stiickelberg fields 0^ as E^^ = 9^i0— 9;^(^ry^j/, when 77^1, = (1,1,1,-1). 
Following a series of arguments^- we can avoiding potential ghost instabilities if we 
fix X^i, = a^ig^^. 

De Sitter solutions are possible, for instance, for ansatz of PG type^ 

ds^ = U^{r,t)[dr + ey/J{r,t)dtf + a^r^{d9^ + sm^ Odip^) - V^{r,t)df. (3) 

In above formula, there are used spherical coordinates labelled in the form u^ = 
[x^ — r,x'^ — 0,y^ = f,y'^ = t), the function / takes non-negative values and 
the constant a = a/{a + 1) and e = ±1. For such bimetric configurations, the 
Stiickelberg fields are parameterized in the unitary gauge a.s (jA — t and (j>^ = 
rn-, (j)^ ~ rn-, (fA — rn-, where a three dimensional (3-d) unit vector is defined 
as n = (n- — sm6 cos (p,n- = sm0sm<f,n- = cos 9). Any PG metric (I3|) defines 
solutions both in GR and in MGT. It allows us to extract the de Sitter solution, in 
the absence of matter, and to obtain standard cosmological equations with FLRW 
metric, for a perfect fiuid source T^^ — [p{t) + p{t)]u^u^ + p(t)g^^, where u^ = 
(0, 0, 0, —V) can be reproduced for the effective cosmological constant '^^f X — fi'^/a. 
It is also possible to express metrics of type ([3]) in a familiar cosmological FLRW 
form (see formulas (23), (24) and (27) in first references'). 
Let us consider an ansatz 

ds^ = 7/1 (r, d)gi ir)dr^ + 772 (r, e)g2{r)dd^ + uj^r, 9, ip, t){r,3{r, 9, t)k{r, 9) (4) 

[d</j + n.,{r, 9)dx^f + r]4{r, 9, t)}n{r, 9, t)[dt + {wi{x'',t) + Wi{x''))dx'f}, 

with Killing symmetry on ^3 — 9^, which (in general) can not be diagonalized 
by coordinate transforms. The values 77^ are called "polarization" functions; uj is 
the so-called "vertical", v, conformal factor. The off-diagonal, N-coefficients, are 
labelled N°-{x'',y*), where (for this ansatz) Nf = ni{r,9) and Nf" = Wi{x'^ ,t) + 
Wi{x^). The data for the "primary" metric are gi[r) = U^ — h4{wi)'^ , g2{r) = 
d?r^M = 5^2 sin^ 6*, /i4 = ^\fU'^ -V\wi = ey/fU^/hi,W2 = 0, n, = 0, when 
the coordinate system is such way fixed that the values /, U, V in ([3]) result in a 
coefficient gi depending only on r. 

We shall work with respect to a class of N-adapted (dual) bases Gq = (e^ = 9,; — 
N^dbjCa = da = d/dy"") and e^ — {e^ — dx'',e'' == dy'' + Nj^dy'^), which are nonholo- 
nomic (equivalently, anholonomic) because, in general, there are satisfied relations of 
type ecG/j — epBa = W^aSj, for certain nontrivial anholonomy coefficients W2aiu). 
For simplicity, we shall consider energy momentum sources and effective cosmologi- 
cal constants which up to frame/coordinate transforms can be parameterized in the 

for constant values ™T and "T with respect to N-adapted frames. 

The system of field equations ^ can be integrated in certain generic off-diagonal 
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general forms depending on all spacetime coordinatesj^l 



ds^ = c'^^^ \{dx^f + {dx^f] + ^ . ^^^ ax) ^3['^'^+(^fc n)dx''f 

-l^YT^^h^dt + (a, I)rfxf . (5) 

for any $ = $, (9,;$)* = di^* and Wi + Wi = 9^$/$* = di A. We can consider 
arbitrary nontrivial sources, ™T+ "T ^ 0, and generating functions, ^{x^,t) := e''^ 
and rife ~ dkn{x'^). The polarization 77-functions for ([5]) are computed in the form 
771 = e^/gi,m = 6^^/52, % = $V4 ( "T + "T), 774 = ($*) V( ™T + "T)$2. 

It is possible to provide an "alternative" treatment of ([5]) as exact solutions 
in MGT. We have to define and analyze the properties of fiducial Stiickelberg 
fields (j)^ and the corresponding bimetric structure resulting in target solutions 
g = {gi,ha,Nj): Let us analyze the primary configurations related to (^^ == ((^ = 
a{T)pQ.~^n-,(jA = a{T)pQ.~^n-,(fA = tk~^), when the corresponding prime PG- 
metric g is taken in FLRW form ds^ = a'^{dp'^ /{l-Kp'^)+p'^{d9'^+sm'^ edip'^))-dT'^. 
The related fiducial tensor is computed Y^^^jdu^u- = ^[dp^ + p^ {dO"^ +s\t? 9dip'^) + 

2H pdpdr — (-^3 — H^ p^)dT^], where the coefficients and coordinates are re-defined 
in the form r ^ p — otr j air) and t ^ t = nt, for K — 0, ±1; k, is an integration 
constant; H := dlna/dr and the local coordinates are parameterized in the form 
X- = p,x- — 0, y- — (p,y- — T. 

In summary, we have found new cosmological off-diagonal solutions in mas- 
sive gravity with flat, open and closed spatial geometries. We applied a geometric 
techniques for decoupling the field equations and constructing exact solutions in 
f{R) gravity, theories with nontrivial torsion and noholonomic constraints to GR 
and possible extensions on (co) tangent Lorentz bundles. A very important prop- 
erty of such generalized classes of solutions is that they depend, in principle, on all 
spacetime coordinates via generating and integration functions and constants. After 
some classes of solutions were constructed in general forms, we can impose at the 
end nonholonomic constraints, cosmological approximations, extract configurations 
with a prescribed spacetime symmetry, consider asymptotic conditions etc. 

Acknowledgments: Authors's research is partially supported by the Program 
IDEI, PN-II-ID-PCE-2011-3-0256. 

References 

1. C. de Rham, Phys. Lett. B 681, 89 (2009); C. de Rham, G. Gabadadze, A. J. ToUey, 
Piys. Rev. Lett. 106, 231101 (2011); S. Hassan, R. A. Rosen, Phys. Rev. Lett. 108, 
041101 (2012); JHEP 1204, 2012 (2012) 

2. P. Stavrinos, S. Vacaru, Class. Quant. Grav. 30, 055012 (2013); S. Vacaru, Gen. Relat. 
Grav. 44, 1015 (2012); EPL 96, 5001 (2011); Class. Quant. Grav. 28, 215991 (2011); 
Int. J. Geom. Meth. Mod. Phys. 8, 9 (2011); Int. J. Mod. Phys. D 21, 1250072 (2012) 

3. T. Kobayashi, S. Siino, M. Yamaguchi and D. Yoshida, Phys. Rev. D 86, 061505 (2012); 
S. Koyama, G. Niz and G. Tasinato, JHEP, 1112, 065 (2011); P. Gratia, W. Hu and 
M. Wyman, Phys. Rev. D 86, 061504 (2012) 



